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ELASTOHYDRODYNAMIC LUBRICATION OF ELLIPTICAL CONTACTS FOR 


MATERIALS OF LOW ELASTIC MODULUS 
I - FULLY FLOODED CONJUNCTION 
by Bernard J. Hamrock and Duncan Dowson* 

Lewis Research Center 

SUMMARY 

Our earlier studies of elastohydrodynamic lubrication of conjunctions of elliptical 
form are applied to the particular and interesting situation exhibited by materials of low 
elastic modulus. By modifying the procedures we outlined in an earlier publication, the 
influence of the ellipticity parameter k and the dimensionless speed U, load W, and 
material G parameters on minimum film thickness for these materials has been inves- 
tigated. The ellipticity parameter was varied from 1 (a ball-on-plate configuration) to 
12 (a configuration approaching a line contact) • The dimensionless speed and load pa- 
rameters were varied by 1 order of magnitude. Seventeen different cases were used to 
generate the following minimum- and central-film-thickness relations: 


H . = 7.43(1 - 0.85 
mm 


H = 7.32(1 - 0.72 e“®*28k)u0.64^0.22 
c 


Contour plots are presented that illustrate in detail the pressure distribution and film 
thickness in the conjunction. 


Professor of Mechanical Engineering, Leeds University, Leeds, England. 



INTRODUCTION 


Only in recent years has the complete numerical solution of the isothermal elasto- 
hydrodynamic lubrication (EHL) of elliptical contacts successfully emerged. The anal- 
ysis requires the simultaneous solution of the elasticity and Reynolds equations. The 
authors' approach to the theoretical solution has been presented in two previous publi- 
cations (refs. 1 and 2), The first of these publications (ref. 1) presents an elasticity 
model in which the conjunction is divided into equal rectangular regions with a uniform 
pressure applied over each region. The second (ref. 2) gives a complete approach to 
the solution of the el as tohydro dynamic lubrication problem for point contacts. 

The most important practical aspect of the EHL point-contact theory (ref. 2) is the 
determination of the minimum film thickness within the contact. That is, the predic- 
tion of a film of adequate thickness is extremely important for the successful operation 
of machine elements in which these thin, continuous, fluid films occur. In a recent 
paper by the authors (ref. 3) the fully flooded results obtained from the theory given in 
references 1 and 2 were presented. A fully flooded condition is said to exist when the 
inlet distance of the conjunction ceases to influence in any significant way the minimum 
film thickness. In reference 3 the influence of the ellipticity parameter and the dimen- 
sionless speed, load, and material parameters on minimum film thickness was inves- 
tigated. Thirty-four different cases were used in obtaining the fully flooded minimum- 
film- thickness formula. 

In reference 4 the basic theory developed in references 1 and 2 was used to study 
the effect of lubricant starvation on the pressure and film thickness within the conjunc- 
tion. A simple expression for the dimensionless inlet boundary distance was obtained. 
This inlet boundary distance defines whether a fully flooded or a starved condition exists 
in the contact. Fifteen cases, in addition to three presented in reference 3, were used 
to obtain simple expressions for the minimum and central film thicknesses in a starved 
conjunction. 

The work presented in references 1 to 4 related to a material of high elastic modu- 
lus (e. g. , steel). The work presented in the present report is for a material of low 
elastic modulus (e.g. , nitrile rubber). For such materials the distortions are large 
even with light loads. Another feature of the EHL of low-elastic-modulus materials is 
the negligible effect of pressure on the viscosity of the lubricating fluid. Engineering 
applications in which elastohydrodynamic lubrication is important for low-eiastic- 
modulus materials include seals, human joints, tires, and elastomeric- material ma- 
chine elements. 

The problem of line contacts, where side leakage of the fluid is ignored, has been 
solved theoretically for low-elastic-modulus materials by Herrebrugh (ref. 5), Dowson 
and Swales (ref. 6), and Baglin and Archard (ref. 7). The solutions of references 5 
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and 6 were obtained numerically and are based on simultaneous solutions of the hydro- 
dynamic and elasticity equations; the analytical solution of reference 7 relied on the 
assumption of a simplified form for the film shape in the contact region. Biswas and 
Snidle (ref. 8) used the approach of reference 7 to solve the point-contact (ball on plate) 
situation. The present work represents, to the best of the authors' knowledge, the first 
attempt at a complete numerical solution of the problem of isothermal elastohydrody- 
namic lubrication of elliptical contacts for low- elastic- modulus materials. In this re- 
port, no assumptions are made as to the pressure or film thickness within the contact, 
and compressibility and viscous effects are considered. 


SYMBOLS 


a 

b 



E 

E' 


F 

G 

H 


0 

H 


c 

min 


I^min 


semimajor axis of contact ellipse 
semiminor axis of contact ellipse 



modulus of elasticity 
2 



normal applied load 

dimensionless material parameter, E'/p^^ 
dimensionless film thickness, h/R^ 

dimensionless central film thickness obtained from EHL elliptical contact 
theory 

dimensionless central film thickness obtained from least- squares fit of data 

dimensionless minimum film thickness obtained from EHL elliptical contact 
theory 

—0 5 

dimensionless film thickness - speed parameter, ‘ 
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dimensionless film thickness - speed parameter from ref. 8 

dimensionless film thickness - speed parameter from present analysis 

dimensionless film thickness - speed parameter from unpublished data by 
W. E. Jamison of Clemson University 

dimensionless film thickness - speed parameter from ref. 9 


Hmln 

Smin 


TT 

min, L 

^^o 

h 

k 

Mp 

P 


Hz 


IV, as 


R 


r 

S 

U 


u 

W 


w 


Hz 


w 


D 


dimensionless minimum film thickness obtained from EHL elliptical contact 
theory for materials of high elastic modulus (ref. 3) 

dimensionless minimum film thickness obtained from least- squares fit of 
data 

dimensionless minimum film thickness for line contact 

constant, initially estimated 

film thickness 

ellipticity parameter, a/b 

0 75 

dimensionless load-speed parameter, WU 
dimensionless pressure, p/E’ 
dimensionless pressure difference 
dimensionless Hertzian pressure 
pressure 

asymptotic isoviscous pressure 

effective radius 

radius of curvature 

natural separation of elliptical solids 

dimensionless speed parameter, ut]q/E'R^ 

mean surface velocity in x-direction, "2 ^ 

o 

dimensionless load parameter, F/E’R 

X 

elastic deformation due to Hertzian pressure distribution 
elastic deformation due to pressure difference, from eq. (2) 
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coordinate systems defined in report 


x, X,xl 

y, Y,Yj 


a. 

V 

^0 

V 

P 

P 

PQ 

Subscripts: 


pressure-viscosity coefficient 
dimensionless viscosity 
atmospheric viscosity 
Poisson's ratio 
lubricant density 
dimensionless density, p/pq 
atmospheric density 


A 

B 


x,y 


solid A 
solid B 

coordinate system defined in report 


THEORETICAL FORMULATION 


The basic theory developed in reference 2 is used here with some minor modifica- 
tions. It was discovered that numerical convergence was considerably better if the di- 
mensionless pressure was written as 


P = P 


Hz 


+ P 


D 


( 1 ) 


%z dimensionless Hertzian pressure 
Pj^ dimensionless pressure difference 

By making use of equation (1), the relevant Reynolds equation can be written as 


j_/^ !!d\ 

9x1 rf ax J 




2 9Y 




9Y 


ax^rf ax j 


+— 9 (pH) 

k2aY\r, dY / \rJ ax 


( 2 ) 


5 



where k = a/b is the ellipticity parameter. Since the Hertzian pressure is known, the 
solution of the Reynolds equation involves solving for Pjj. In equation (2), compressi- 
bility and viscous effects are considered even though it was mentioned in the 
INTRODUCTION that the effect of pressure on viscosity and density is negligible. 

The equation for the dimensionless film thickness is written as 




+ 


S + w 


Hz 



R 


X 


( 3 ) 


constant (initially guessed) 
natural separation of elliptical solids 

elastic deformation due to Hertzian pressure distributions 
elastic deformation due to pressure difference obtained from eq. (2) 

The evaluation of the elastic deformation is exactly (hat used in references 1 and 2. 

The nodal structure used for obtaining all the results is shown in figure 1. It 
should be noted from this figure (hat, because of the dimensionless representation of the 
coordinates, the actual Hertzian contact ellipse becomes a circle regardless of the value 
of k. The nodal structure shown in figure 1 was arrived at after much exploration in 
which the number of nodes in the semimajor and semiminor axes, as well as the dis- 
tance from the center of the contact to the edges of the computing zone, was varied. 


H, 


w- 


Hz 


w 


D 


DIMENSIONL-ESS GROUPING 

The same dimensionless grouping of parameters used in the analysis of elastohy- 
drodynamic lubrication of high- elastic- modulus materials (ref. 3) is used herein: 

Dimensionless film thickness: 


E=-^ (4) 

I'x 

Ellipticity parameter: 

k=:2 (5) 

b 
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Dimensionless speed parameter: 


U = 




Dimensionless load parameter; 


W = 


F 


E'r| 


Dimensionless material parameter: 


G =■ 


E» 


IV, as 


( 6 ) 


( 7 ) 


( 8 ) 


where p.^ is the asymptotic isoviscous pressure obtained from Roelands (ref. 10). 
The asymptotic isoviscous pressure can be approximated by the inverse of the pressure- 
viscosity coefficient, p-^ 

The dimensionless film thickness can be written as 


H = f(k,U,W, G) 


( 9 ) 


The most important practical aspect of the elastohydrodynamic lubrication elliptical 
contact theory is the determination of the minimum film thickness within the conjunc- 
tion. Therefore, in the fully flooded results presented herein, the dimensionless pa- 
rameters k, U, W, and G were varied, and the effect upon minimum film thickness 
studied. Care was taken to ensure that all results were in the elastic region. 


EFFECT OF ELLIPTICITY OF ELASTIC CONJUNCTION 

The ellipticity parameter is a function of the radii of curvature of the solids only 
(r^, rgj^, r^y, and r^y). The radii of curvature in the x-direction for both solids A 
and B are used in defining the dimensionless speed and load parameters. Therefore, 
only the radius of curvature of solid B in the y-direction was changed in varying k 
from 1 (a ball-on-plate configuration) to 12 (a configuration approaching a line contact). 
In doing this, the dimensionless speed, load, and material parameters were held con- 
stant at 
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U = 0.1028x10' 


-7 




,-3 


W = 0. 4405x10 > 


G = 0.4276 


( 10 ) 


Table I gives seven values of the ellipticity parameter and the corresponding mini- 
mum film thickness as obtained from the EHL elliptical contact theory. Having these 
seven pairs of data, the object was to determine an equation that would describe how the 
ellipticity parameter affects the minimum film thickness. The general form of this 
equation can be written as 



H 


min 


H 


min, L , 


Ae 


Bk 


A leasts squares exponential curve fit to the seven pairs of data points 


ki, 


H 


1 - 


mm 


H 


min, L/ . 


i = 1 , ..., 7 


( 11 ) 


was used in obtaining values for A and B in equation (11). Besides a least-squares 

fit, a coefficient of determination r was obtained. The value of r reflects the fit of 

the data to the resulting equation: 1 being a perfect fit, and zero the worst possible fit. 

The minimum film thickness for a line contact H„. ^ used in equation (11) was de- 

niiii, ij 

termined bv finding the H . t that gives a coefficient of determination closest to 1. 

. mm, L _0 

This value of H . t turned out to be 240. 0x10“ , with a corresponding coefficient of 

mill, Ij 

determination of 0. 9908, which is an excellent fit. Furthermore, the values of A 
and B in equation (11), as obtained from the least- squares fit, are 


A = 0.8488 « 0.85 

( 12 ) 

B = -0.3115 » -0.31 

(13) 


From equations (11) to (13), the following equation can be written, which shows the effec 
effect of k on 0 
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tnin 


/-. n oc -0.31k, 
tx (1 - 0. 85 e ) 


(14) 


Jj’lgure 2 shows the variation of the ratio /H . _ with k for the EHL high- 

el asticj- modulus analysis (ref. 3) and the EHL low-elastio-modulus analysis (present 
results). If we assume that the minimum film thickness obtained from the EHL ellipti- 
cal analysis can only be obtained to an accuracy of 3 percent, we find a limiting solution 
for line contact being approached for a k of 5 for the EHL high- elastic- modulus analy- 
sis and for a k of 11 for the EHL low-elastic-modulus analysis. 


INFLUENCE OF SPEED 


Changing only the surface velocity in the x- direction u causes the dimensionless 
speed parameter U to change while the other dimensionless parameters (k, W, and G) 
remain constant. The values at which these dimensionless parameters were held con- 
stant in the calculations were 


W = 0.4405x10"^ 

k = 6 r 


(15) 


G = 0.4276 


Table I also gives the dimensionless speed parameter and the corresponding di- 
mensionless minimum film thickness as obtained from the EHL elliptical contact theory. 
There are five different values of the dimensionless speed parameter covering 1 order 
of magnitude. Having these five pairs of data, the objective was to determine an equa- 
tion that would describe how the speed affects the minimum film thickness. The general 
form of this equation can be written as 


H • 
min 


= lU 


J 


(16) 


By applying a least-squares power fit to the five pairs of data, (U,-,H ^.^ .), where 
i = 1, . . . , 5, the values of I and J were found to be 


1= 32.48 


(17) 


9 



J = 0.6505 w 0.65 


(18) 
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The coefficient of determination r for these results was excellent at 0. 9997. From 
equations (16) and (18) the effect of dimensionless speed on dimensionless minimum 
film thickness can be written as 


H . oc U 
mm 


Figure 3 gives contour plots of dimensioxiless pressure for two extreme values of 
—7 —7 

U, 0.05139x10 and 0. 5139x10 . In these and all contour plots to be presented, the 

+ symbol indicates the center of the Hertzian contact. Because of the dimensionless 
representation of the X and Y coordinates, the actual Hertzian contact ellipse be- 
comes a circle regardless of the value of k. The Hertzian contact circle is shown in 
each figure by asterisks. At the top of each figure the contour labels and the corre- 
sponding values of dimensionless pressure are given. The inlet region is to the left, 
and the exit region to the ri^t. 

The pressure contours shown in figure 3 are nearly circular, or Hertzian. In fig- 
ure 3(b), the high-speed case, the inlet pressure was hi^er than in the low-speed case 
shown in figure 3(a). Inside the contact the contour values of the dimensionless pres- 
sure were higher for the low- speed case (fig. 3(a)). The pressure spikes found when 
dealing with materials of high elastic modulus (ref. 3) were not in evidence in these so- 
lutions for low - elastic- modulus materials. The lack of a pressure spike is due to the 
lack of viscous effects of the fluid in the contact for low-elastic-modulus materials. 
This in turn is due to the fact that considerably less pressure is generated in a contact 
with low-elastic-modulus materials than in a contact with high- elastic- modulus ma- 
terials. 

Figure 4 shows contour plots of dimensionless film thickness for two values of 

-7 -7 

U, 0.05139x10 and 0.5139x10 . Figure 4(a) shows two regions of minimum film 

thickness, close to the Hertzian circle and off to the side. Figure 4(b) shows one 
region of minimum film thickness, between the center of the contact and the Hertzian 
circle. 

Figure 5 shows the variation of pressure and film thickness on the X-axis near the 
midplane of the conjunction for the same two values of U. For all the solutions for 
various speeds, the values of the dimensionless load, material, and ellipticity parame- 
ters were held fixed as described in equation (15). In figure 5 the pressure in the inlet 

“■ 7 

region is higher for the higher speed (U = 0. 5139x10 ) profile. It should be noted that 

these sli^t changes in pressure that occur when the speed is increased by 1 order of 
magnitude result in a significant change in film thickness, as indicated by equation (19). 


10 



This illustrates most clearly the dominant effect of the dimensionless speed parameter 
on the minimum film thickness in elastohydrodynamic contacts for low-elastic-modulus 
materials. Similar results were found for hi^- elastic- modulus materials (ref. 3). 


INFLUENCE OF LOAD 


Changing only the normal applied load F in equation (7) causes the dimensionless 
load parameter W to change while the other dimensionless parameters (k, U, and G) 
remain constant. The values at which these parameters were held constant were 


U = 0.1028x10' 


,-7 


k= 6 


G = 0.4276 


( 20 ) 


The load results presented in table I give the dimensionless load parameter and the 
corresponding minimum film thickness as obtained from the EHL elliptical contact 
theory. There are six different values of dimensionless load parameter covering 1 
order of magnitude. Having these six pairs of data, the objective was to determine an 
equation that would describe how the dimensionless load affects the minimum film thick- 
ness. The general form of this equation can be written as 


C = KW^ 


( 21 ) 


where 


C = 


Hmin 

(1 - 0.85 


( 22 ) 


In equation (22) the exponents are rounded off to two significant figures so that any error 
could be absorbed in K, given in equation (21). By applying a least- squares power fit 
to the six pairs of data [ (W., C.), i = 1, . . . , 6], the values of K and L were found 
to be 


K = 6.507 


(23) 
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L = -0.2075 « -0.21 


(24) 
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The coefficient of determination r for these results was 0. 9985, which is excellent. 
From equations (21), (22), and (24), the effect of load on minimum film thickness can 
be written as 


mm 


oc W“ 


0.21 


(25) 


Figure 6 gives contour plots of dimensionless pressure for the two extreme values 

-3 -3 

of W that were investigated, 0.2202x10 and 2.202x10 . Again the pressure con- 

tours were nearly circular, or Hertzian. The contour values were considerably higher 
for the high-load case (fig. 6(b)). 

Contour plots of dimensionless film thickness for the same two values of W are 

—3 

shown in figure 7. In figure 7(a), for the low-load case (W = 0.2202x10 ), the mini- 

mum film thickness occurred directly behind the center of the contact; in figure 7(b), 

-3 

for the hi^-load case (W = 2.202x10 ), the minimum film thickness also occurred di- 

rectly behind the center of the contact but closer to the Hertzian circle. The two C 
contours in figure 7(b) indicate a sli^t increase in film thickness before the minimum- 
film- thickness region is reached. 

The variation of pressure and film thickness in the X-direction along a line close 
to the midplane of the conjunction is shown in figure 8 for two values of W. The 
values of U, G, and k were held fixed as described by equation (20) for all computa- 
tions at various loads. In this figure, for a 1- order- of- magnitude change in the dimen- 
sionless load parameter, there is a considerable change in the pressure profile but not 
such a significant change in the film thickness profile. 


MINIMUM- FILM- THICKNESS FORMULA 

The proportionality expressions (14), (19), and (25) established how the minimum 
film thickness varied with the ellipticity, speed, and load parameters, respectively. 
This enabled a composite minimum-film-thickness formula for a fully flooded, iso- 
thermal, elastohydrodynamic elliptical contact for low-elastic modulus materials to be 
written as 


H . = 7.43(1 - 0.85 (26) 
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In equation (26) the constant 7.43 is different from that given in equation (23) to account 
for rounding off the load-parameter exponent. 

Table I gives the values for the minimum film thickness obtained from the least- 
squares fit as defined by equation (26). The percentage difference between the minimum 
film thickness obtained from the EHL elliptical contact theory and the minimum 

film thickness obtained from the least- squares fit equation H , is expressed as 


D 


1 


/H . - H . , 
_/ mm mm 


100 


H 


mm 


(27) 


In table I the values of are within the range -8 and +3. 

It is sometimes more convenient to express the side-leakage factor in equation (26) 
in terms of the curvature ratio Ry/Rx ^stead of the ellipticity parameter through the 
following relation: 


k = 


1.03 



where 




’^Ay ^By 



(28) 


(29) 


(30) 


Using equation (28) avoids the need to evaluate elliptic integrals of the first and second 
kinds in the determination of k. The minimum film thickness can thus be derived di- 
rectly from a knowledge of the curvature of the contacting bodies (r^, r^^, r^^, 
and r^y). 

It is interesting to compare the equation for materials of low elastic modulus 
(eq. (26)) with the corresponding equation generated in reference 3 for materials of hi^ 
elastic modulus 


mm 


(31) 
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The powers of U in equations (26) and (31) are quite similar, but the power of W is 
much more significant for low- elastic- modulus materials. The expression showing the 
effect of the ellipticity parameter is of exponential form in both equations, but with quite 
different constants. 

A major difference in equations (26) and (31) is the absence of a material parameter 
in the expression for the minimum film thickness for low- elastic- modulus materials. 
There are two reasons for this. One is the negligible effect of pressure on the viscos- 
ity of the lubricating fluid for low-elastic-modulus materials. The other is the way in 
which the role of elasticity is simply and automatically incorporated in the prediction 
of conjunction behavior through an increase in the size of the Hertzian contact zone cor- 
responding to changes in load. As a check on the validity of this, case 9 of table I was 
repeated with the material changed from nitrile rubber to silicone rubber. This is re- 
corded as case 17 (table I). As can be seen from equations (6) to (8), when the ma- 
terial (as expressed by E') is changed, not only does the material parameter change, 
but so do the dimensionless speed and load parameters. Only the ellipticity parameter 
can be held fixed. From table I, case 17, we find that the minimum film thickness from 
the EHL elliptical contact theory was 181.8x10" . The minimum film thickness from 
the least-squares fit (eq. (26)) turned out to be 182. 5x10“ . This clearly indicates a 
lack of dependence of the minimum film thickness for low-elastic- modulus materials on 
the material parameter. 


CENTRAL- FILM- THICKNESS FORMULA 

There is interest in knowing the central film thickness, in addition to the minimum 
film thickness, in elastohydrodynamic contacts. The procedure used in obtaining the 
central film thickness was the same as that used in obtaining the minimum film thick- 
ness and is not repeated here. The central-film-thickness formula for low-elastic- 
modulus materials obtained from the results is 

H^ = 7.32(1 - 0.72 (32) 


Comparison of the central-film- thickness formula (eq. (32)) with the minimum-film- 
thickness formula (eq. (26)) reveals only slight differences. 

Table II gives the 16 different cases used to obtain equation (32), as well as case 17, 
which is a check on equation (32). In this table, H^ corresponds to the central film 
thickness obtained from the EHL elliptical contact theory and H^ corresponds to the 
central film thickness obtained from equation (32). The percentage difference between 
these two values is given by Dg and is written as 
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( 33 ) 



In table II the values of Dg are within the range -11 to +23. 

The ratio of minimum to central film thickness evident in the computed values 
ranged from 70 to 83 percent, the average being 77 percent. 


Comparison of Different Investigators* Results 

To evaluate the dimensionless minimum-film- thickness equation (eq. (26)) de- 
veloped in this report, a comparison was made between it, the numerical solution ob- 
tained by Biswas and Snidle (ref. 8), and the recent experimental findings of Jamison 
(unpublished data). Both are applicable only for k = 1. 

From reference 8 the dimensionless minimum film thickness can be written as 


where 


Hmin = «min ^ 


-0.5 


(35) 


Mp = WU 


-0.75 


(36) 


The dimensionless groups given in equations (35) and (36) were first used by Moes and 
Bosma (ref. 11). 

The central-film- thickness equation obtained from the experimental results of 
Jamison can be written as 




= 2.4 


0.075 


(37) 


where 


U 

c c 


-0.5 


(38) 
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I 


mill II 


Modifying equation (37) to provide a minimum-film-thickness equation by using the as- 

yields 

Equations (34) and (39) are valid only for k = 1. 

By making use of equations (35), (36), and table I equation (26) can be written as 

^min} = 8.53(1 - 0.85 (40) 

Therefore, for k= 1, equation (40) reduces to 


sumption that {Hmin}j " {^c}j 


H 


mm 


kFl 1 


= 3.21 M' 


- 0.21 


(41) 


H 


Note the larger negative exponent in equation (41) than in either equations (34) or (39) . 

Figure 9 compares the different investigators’ results from equations (34), (39), 
and (41) . The three equations seem to agree quite well with each dther. The present 
report's results are equivalent to the Biswas and Snidle theoretical results (ref. 8) at 
Mp = 100 and the Jamison experimental results at Mp = 45. Therefore, even though 
the exponent on Mp for the present results is larger than those obtained in reference 8 
and by Jamison, the agreement is quite good for k = 1. Also shown in figure 9 is the 
rigid isoviscous solution obtained from Kapitza (ref. 9). 

Figure 10 uses equation (40) to show the effect of the dimensionless parameter Mp 
on the dimensionless minimum film thickness for six values of the ellipticity parameter. 
As k increases, less change occurs in the dimensionless minimum film thickness. 
From this figure the dimensionless minimum film thickness can easily be obtained for 
known values of k and Mp. 

Similarly, the equation for the dimensionless central film thickness can be written 
as 


|h^| = 11.61(1 - 0.72 


(42) 
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Figure 11 uses equation (42) to show the effect of the dimensionless parameter Mp on 
the dimensionless central film thickness for six values of k. 


CONCLUDING REMARKS 


By modifying the procedures outlined by the authors in an earlier publication, the 
influence of the ellipticity parameter k and the dimensionless speed U, load W, gnri 
material G parameters on minimum film thickness for low- elastic- modulus materials 
has been investigated. The ellipticity parameter was varied from 1 (a ball-on-plate 
configuration) to 12 (a configuration approaching a line contact) . The dimensionless 
speed and load parameters were varied by 1 order of magnitude. Seventeen different 
cases were used to generate the minimum- and central-film-thickness relations: 


Hmin= 7.43(1 - 0.85 


= 7.32(1 - 0.72 


g-0.28k^^j0.64^0.22 


Contour plots have been presented that show in detail the pressure distribution and the 
film thickness in the conjunction. 

The present report presents for the first time a complete theoretical film- thickness 
solution for the problem of isothermal elastohydrodynamic elliptical contacts for low- 
elastic- modulus materials operating under fully flooded conditions. 
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TABLE I. - DATA SHOWING EFFECT OF ELLIPTICITY, LOAD, SPEED, AND MATERIAL PARAMETERS ON MINIMUM FILM THICKNESS 


Case 

Ellipticity 

parameter, 

k 

Dimensionless 

load 

parameter, 

W 

Dimensionless 

speed 

parameter, 

U 

Dimensionless 

material 

parameter, 

G 

Minimum film thickness 

Difference 

between 

Hmin H^i„. 

Di. 

percent 

Results 

Obtained from 
EHL elliptical 
contact theory, 

min 

Obtained from 
leasts squares 
fit, 

min 

1 


1 


0.1028x10“'^ 

0.4276 

88.51x10"® 

91.08>d0"® 

+2.90 



2 


2 








142.5 

131.2 

-7.93 



3 


3 








170.4 

160.8 

-5.63 



4 


4 








186.7 

182.4 

-2.30 


r Ellipticity 

5 


6 








206.2 

209.8 

+1.75 



6 


8 








219.7 

224.6 

+2.23 



7 

12 








235.2 

236.0 

+. 34 



8 


6 



.05139 



131.8 

133.7 

+1.44 



9 





.1542 



268.1 

273.1 

+1.86 

















^Speed plus case 5 






.2570 



381.6 

380.7 

-.24 



11 





.5139 



584.7 

597.3 

+2.15 

J 


12 



.2202 

.1028 



241.8 

242.7 

+. 37 



13 



.6607 






190.7 

192.7 

+1.05 



14 










170.5 

173.1 

+1.52 


^Load plus case 5 

15 



1.542 







160.4 

161.3 

+. 56 



16 





T 




149.8 

149.7 

-.07 



17 



.1762 

.06169 

1.069 

181.8 

182.5 

+.39 


Material 












to 
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TABLE II. - DATA SHOWING EFFECT OF ELLIPTICITY, LOAD, SPEED, AND MATERIAL PARAMETERS ON CENTRAL FILM THICKNESS 



Ellipticity 

parameter, 

k 


Dimensionless Dimensionless Dimensionless 


Minimum film thickness 


load 

parameter, 

W 


0.4405x10" 


speed 

parameter, 

U 


0.1028x10" 


8 

6 


.05139 

9 




.1542 

10 




.2570 

11 



If 

.5139 

12 



.2202 

,1028 

13 



.6607 


14 



1.101 


15 



1.542 


16 



2.202 


. 

1 

if 

.1762 

i 

.06169 


material 

parameter, 

G 


Obtained from 
EHL elliptical 
, contact theory. 


Obtained from 
least squares 


Difference 

between 

H^ and H^ 
c c 


114. 9x10"® 

141.0x10® 

+22.72 

204.2 

182.2 

-10.77 

217.2 

213.2 

-1.84 

253.1 

236.7 

-6.48 

265.0 

267.9 

+1.09 

288.8 

285.7 

-1.07 

307.1 

301.7 

-1.76 

171.7 

171.9 

+. 12 

354.0 

347.3 

-1.89 

498.7 

! 481.5 

-3.45 

743.7 

! 750.3 

+.89 

316.3 

312.0 

-1.36 

235.1 

245.0 

+4.21 

221.5 

219.0 

-1.13 

192.6 

203.4 

+5.61 

196.6 

188.0 

-4.37 

237.4 

236.3 

.46 



Ellipticity 


Speed plus case 5 


Load plus case 5 


Material 







Figure 1. - Nodal structure used for numerical calculations. 


Film thickness ratio. 



0 2 4 6 8 10 

Ellipticity parameter, k 


Figure 2, - Variation of ratio of dimensionless minimum film thick- 
ness to dimensionless minimum film thickness for a line contact 
with ellipticity parameter, for EHL high- and low-elastic-modulus 


to 

to 


Dimensionless pressure, 
P =p/E' 

A 1.285x10‘2 

B 1.270 

C 1.220 

D 1.120 

E .970 

F .770 

G .520 

H .200 



Dimensionless pressure, 
P =p/E' 

A 1.230x10'2 
B 1.210 

C 1.170 

D 1.100 

E 1.000 

F .850 

G .600 

H .300 



(a) Contour plot of dimensionless pressure for dimensionless speed U of (b) Contour plot of dimensionless pressure for dimensionless speed U of 

a 05139xl0’^. 0. 5139x10"’ 

Figure 3. - Contour plots of dimensionless pressure for dimensionless speed parameters U of 0. 05139x10"^ and 0. 5139x10"^. The dimensionless parameters 
k, W, and G are held constant as defined in equation (15). 



Dimensionless 
film thickness, 
H«h/Rx 
A 1.35x10*’^ 

B 1.40 
C 1.50 
D 1.60 
E 1.80 
F 2.10 
G 2.50 
H 3,00 



(a) Contour piot of dimensionless film thickness for dimensionless speed U 
of 0,05139x10-7. 


Dimensionless 
film thickness, 
H » h/Rx 
A 5.88x10-4 
B 6.00 
C 6.20 
D 6.40 
E 6.70 
F 7.20 
G 8.00 
H 9.20 



Figure 4. - Contour plots of dimensionless film thickness for dimensionless speed parameters U of 0, 05139xl0"7 and a 5139x10-7. The dimensionless 
parameters k, W, and G are held constant as defined in equation (15). 


Dimensionless film thickness, H « h/R 











Number of nodes in X-direction 


Figure 8. - Variation of dimensionless pressure and film thickness on X-axis for dimensionless load parameters W of 
0. 2202x10’^ and 2.202x10"^. The value of Y is held fixed near axial center of contact. 



Dimensionless load-speed parameter. Mp = 

Figure 9. - Comparison of different investigators* results. 
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Figure 10. - Effect of dimensionless load-speed parameter on dimensionless minimum 
film thickness - speKi parameter for six eliipticify parameter values. 



Figure 11. - Effect of dimensionless load-s^ed parameter on dimensionless central film 
thickness - speed parameter for six ellipticify parameter values. 


28 


NASA-Langley, 1977 


E-9216 


Illlllllllllll I II 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTV for PRJVATE use *300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AtsTD FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



bio 001 Cl 0 D 770715 S00903DS 
DEPT OF THE 8IB FCHCE 
AF WEAPONS LIEOBATORI 
ATTN: TECHNICAL LIBRARY (SDL) 
KIRTIAND AFE NH 87117 


POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


aeronautical and space activities of the United States shall be 
conducted so as to contribute , . . to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


•neCHNICAL REPORTS: Sdencific and 

technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obfained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



